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Two-particle dark state cooling of a nanomechanical resonator 
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The steady-state cooling of a nanomechanical resonator interacting with three coupled quantum dots is stud- 
ied. General conditions for the cooling to the ground state with single and two-electron dark states are obtained. 
The results show that in the case of the interaction of the resonator with a single-electron dark state, no cooling 
of the resonator occurs unless the quantum dots are not identical. The steady-state cooling is possible only 
if the energy state of the quantum dot coupled to the drain electrode is detuned from the energy states of the 
dots coupled to the electron source electrode. The detuning has the effect of unequal shifting of the effective 
dressed states of the system that the cooling and heating processes occur at different frequencies. For the case 
of two electrons injected to the quantum dot system, the creation of a two-particle dark state is established to 
be possible with spin-antiparallel electrons. The results predict that with the two-particle dark state, an effective 
cooling can be achieved even with identical quantum dots subject of an asymmetry only in the charging potential 
energies coupling the injected electrons. It is found that similar to the case of the single-electron dark state, the 
asymmetries result in the cooling and heating processes to occur at different frequencies. However, an important 
difference between the single and two-particle dark state cases is that the cooling process occurs at significantly 
different frequencies. This indicates that the frequency at which the resonator could be cooled to its ground state 
can be changed by switching from the one-electron to the two-electron Coulomb blockade process. 

PACS numbers: 42.50.Wk, 42.50.Lc, 07.10.Cm 



I. INTRODUCTION 



Cooling of a nanomechanical resonator (NMR) to its quan- 
tum ground state is a field of interesting theoretical and exper- 
imental research owing to their potential applications in high- 
precision detection of mass |[T|], mechanical displacement 121, 
and quantum information processing iJ-lsl]. One of the con- 
ventional methods is the sideband cooling approach, based on 
which various cooling schemes have been put forward by cou- 
pling the resonator to a dissipative two-level system such as 
superconducting qubits |6, 7], nitrogen- vacancy defects in di- 
amond 1^, quantum dots ^ [lOll and suspended carbon nan- 
otubes lllll[l2ll . Another method that could lead to an efficient 
cooling to the ground state involves the destructive quantum 
interference, i.e. the creation of a dark state that suppresses 
the carrier and blue sideband transitions. The method has 
inspired numerous cooling schemes for trapped ions such as 
electromagnetically induced transparency (BIT) cooling I113fl 
and Stark-shift cooling 1141] . Besides these, combining the 
EIT cooling scheme with Stark-shift cooling scheme, a robust 
and fast laser cooling for trapped ions, atoms or cantilevers 
with closed interaction contour has been predicted [ISl. It has 
also been demonstrated that cooling of a single NMR can be 
achieved via dynamical back action process or radiation pres- 
sure acting on the NMR lfl6l[l7ll . 

There have recently been a number of theoretical studies of 
the transport properties of a triple quantum dot (TQD) sys- 
tem arranged in a triangular geometry lIlSl - EOll . It has been 
demonstrated that by tuning of the ac field frequency, the anti- 
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resonant behaviour in the current occurs. As the result, the 
stationary current is found to be sensitive to A-B effect due 
to spin-blockade [18]. In addition, the role of the dark state 
has been considered with two of the dots weakly connected 
to external leads lll9ll . The results showed that the coupling 
between the dot electrons and the dissipative single-phonon 
mode can influence the transport properties of the system [20|]. 
However, it has been found that the presence of the dark state 
could be important not only in the study of transport proper- 
ties but also for the cooling of micro mechanical or nanome- 
chanical resonators. In particular, Li et al. 12111 have recently 
demonstrated that an NMR can be cooled to its ground state by 
coupling to the TQD system operating in the strong Coulomb- 
blockade regime [1221 • 

As we have already mentioned, the cooling schemes con- 
sidered so far have been limited to the case of the strong 
Coulomb-blockade regime, in which at most one electron is 
allowed in the cooling system at a given time. However, 
the strong Coulomb-blockade situation may not always work 
since a slight change in the chemical potential of the elec- 
trodes could lead to the inclusion of two electrons into the 
transport window. With the possibility to inject two electrons, 
one can expect to see novel features in the electron transport 
and cooling properties of a given system. For example, Poltl 
et al. 1I23I1 have shown that the presence of two electrons in 
the TQD system may lead to the creation of a two-electron 
dark state, which could result in breaking of the transport 
current and the system then exhibiting super-Poissonian be- 
haviours 1I23I1 . 

It is the purpose of this paper to study the cooling proper- 
ties of the TQD system when two electrons could be present in 
the dot system. Comparison is made with the single-electron 
cooling calculated by Li et al. UM that arises in the strong 
Coulomb-blockade regime. We include charging potential en- 



ergies to describe the inclusion of two electrons. Following 
the prediction of Poltl et al. 12311 that under some circum- 
stances a two-electron dark state can be created in the dot 
system, we demonstrate how the presence of this dark state 
could lead to cooling of the NMR to its ground state. We find 
that in both, the single and two-electron dark state cases, the 
TQD system may lead to cooling of the NMR only if there is 
an asymmetry present in the system. In particular, in the case 
of single electron injected to the dot system, a cooling of the 
NMR occurs only when the energy state of the quantum dot 
couple to the drain electrode is detuned from the energy states 
of the dots coupled to the electron source electrode. On the 
other hand, with two electrons present in the dot system, cool- 
ing of the NMR can be achieved even with identical quantum 
dots subject of an asymmetry between the charging potential 
energies coupling the injected electrons. 

The paper is organized as follows. In Sec.|II]we introduce 
the setup and explain the Hamiltonian of the whole system, 
then we derive the master equation and the rate equations for 
the cooling dynamics. A detailed calculation of the conditions 
for cooling of an NMR with single- and two-electron dark 
states is presented in Sec. |lll] We summarize our results in 
Sec. HV] Finally, transformation matrices from the bare elec- 
tron states to the dressed electron states are given in the Ap- 
pendixlAl 



II. DESCRIPTION OF THE SYSTEM 

A schematic diagram of the system we consider is shown 
in Fig. [T] Three quantum dots, each composed of a single- 
electron orbital state, are arranged in a triangular geometry 
and are weakly connected to three fermionic electrodes via 
tunnel barriers. Two dots, labeled 1 and 2, are capacitively 
coupled to each other without electrons tunnelling directly be- 
tween them. The dots 1 and 2 are coupled to the dot 3 through 
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FIG. 1 : Schematic diagram of the system composed of three quantum 
dots connected to the source and drain electrodes through tunneOing 
baniers. The dots 1 and 2 are capacitively coupled to each other 
without electrons tunnelling directly between them while both dots 1 
and 2 are coupled coherently to the dot 3 with the tunnelling ampli- 
tudes Ti and T2, respectively. A nanomechanical resonator (NMR) 
is capacitively coupled to the dots 2 and 3 only. 



the electron tunnelling effect that occurs with amplitudes Ti 
and T2, respectively. Thus, electrons can tunnel from the dot 1 
to 3 and from 2 to 3 with no tunnelling permitted from the 
dots 1 and 2. The cross correlation between two capacitively 
coupled quantum dots in the Coulomb blockade regime has 
recently been measured ll24ll . We assume that only two quan- 
tum dots, 2 and 3 are directly coupled to the NMR. 



A. Hamiltonian of the system 

We can divide the Hamiltonian of the entire system into 
four parts: (1) Three coupled quantum dots, (2) three normal 
metal electrodes coupled to the quantum dots, (3) a mechan- 
ical resonator interacting with two of the three quantum dots, 
and (4) a reservoir or thermal bath surrounding the mechanical 
resonator 

The dynamics of the three coupled quantum dots contain- 
ing maximum two electrons can be determined by the Hamil- 
tonian 



i=l a a 
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(1) 



where nia = dl^dia is a number operator describing the num- 
ber of spin-up {a =t) ^nd spin-down (cr =4,) electrons in 
the ith quantum dot. The operators di^ and d]^ are respec- 
tively the annihilation and creation operators of an electron 
with spin a in the ith quantum dot, and Ei is the energy of 
the single-electron state in the dot i. We assume that the 
energy states of the electron in each quantum dot are spin- 
independent, so both spin-up and spin-down electrons can 
equally occupy the energy states of the quantum dots. The 
coupling between the quantum dots is arranged through a tun- 
nelling process such that the tunnelling of electrons between 
the dots 1 and 3 occurs with an amplitude Ti, whereas be- 
tween the dots 2 and 3 occurs with an amplitude T2 . There is 
no direct tunnelling allowed between the dots 1 and 2. In this 
case, the system may be considered as a Lambda-type system 
with two ground states 1 and 2 and the upper state 3. 

A possibility of simultaneous injection of two electrons into 
the dots is determined in the Hamiltonian ([T]i by charging po- 
tential energies. There are two distinct types of the charging 
potential energies Uu and Uij [i ^ j), corresponding to the 
charging energy of two electrons injected to the same ith dot 
and to different dots, respectively. Notice that the charging 
potential energies Uij are independent of the polarization of 
the spins. 

The electrons can be injected to the quantum dots from nor- 
mal metal electrodes by coupling them through a tunnelling 
process. Suppose that the electrodes coupled to the dots 1 
and 2 appear as a source (donors) of the electrons while the 
electrode coupled to the dot 3 appears as an absorber (drain) 
of the electrons. In this case, the Hamiltonian of the electrodes 



coupled to the quantum dots can be written as 



B. One and two electron dynamics 






Cika + [Vikcl.d 



jk^jka"']'^ 



H.c. 



, (2) 
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where the first term is the energy of the electrons with ejk 
representing the energy of a single electron of wave number k 
in the electrode j. The operators Ci^a and c\j^^ are respec- 
tively the annihilation and creation operators of an electron 
with spin a and wave number k in the electrode j. The sec- 
ond term in Eq. (|2]l is the interaction between the electrodes 
and the quantum dots with Vjk denoting the coupling strength 
between an electron of the wave number k in the jth electrode 
and the corresponding quantum dot. The coupling is arranged 
through a tunnelling of the electrons between the electrodes 
and the quantum dots, and is independent of a, so that both 
spin-up and spin-down electrons are coupled equally to the 
quantum dots. 

The capacitor adjacent to the dots is formed by the NMR 
and a static plate with a gate voltage applied to it. The effec- 
tive capacitance of this capacitor could be adjusted by the dis- 
placement of the NMR from its equilibrium position lUOl 12111 . 
In this way, the displacement of the NMR couples to the num- 
ber of electrons in the dots. We assume that only two of them, 
the dots 2 and 3 participate in the interaction with the NMR. 
In this case, the Hamiltonian of the NMR interacting with the 
quantum dots is of the form Ill2ll25ll 

Hyn = ujma)a + V] (a + o)) {ctifiia + a2"-3o-) , (3) 



where a^ and a are operators representing the creation and 
annihilation of a phonon excitation of frequency uj„i in 
the NMR, and a^ (i = 1, 2) is the electromechanical coupling 
constant between the ith quantum dot and the phonon mode. 
Note that the interaction of the NMR with the dots is in a form 
of the nonlinear parametric interaction, known in optical sys- 
tems as a radiation pressure type interaction [26]. 

In addition to the interaction with the quantum dots, the 
phonon mode of the NMR is also allowed to interact with a 
heat bath or thermal reservoir which results in the presence of 
thermal phonons in the NMR mode. The Hamiltonian for this 
case is simply written as 



m 



= E' 



^t; 



^9q{ala + aqa) 



(4) 



where al and dq are respectively the creation and annihilation 
operators of the qth mode of frequency Uq of the multi-mode 
reservoir, and gq is the coupling strength of the mode q to 
the NMR phonon mode d. 

We should point out here that in the Hamiltonians ([Tli-(|4]i 
the coupling parts retain only the terms which play a dom- 
inant role in the rotating-wave approximation. Antiresonant 
terms which would make much smaller contributions have 
been omitted. 



The possibility of a simultaneous injection of two electrons 
to the quantum dots system can be directly addressed by writ- 
ing the creation operators d\^ in the two-electron basis as 



< = 


|lt0)(00| + |lt2t>(2t0| + 


|lt3t)(3t0| + 


ltli)(UO| 


+ 


|lt2^)(2^0| + |lt3^)(3^0|, 






4t = 


|2^0)(00|-|lt2t)(lt0| + 


|2t3t)(3t0| + 


2tU)(U0| 


+ 


|2t2j(2^0| + |2t3;)(3;0|, 






4t = 


|3t0)(00|-|lt3t)(lt0|- 


|2t3t)(2t0| + 


3TU)(liO| 


+ 


|3t2;)(2;0| + |3t3;)(3;0|, 






d\, = 


|UO)(00|-|ltU>(ltO|- 


|2tU>(2tO|- 


3tU)(3tO| 


+ 


|lA>(2;0| + |lA>(3iO|, 






dl- 


|2^0)(00|- 11^2^(1^01- 


|2t2;)(2t0|- 


3T2i>(3tO| 


- 


|1A>(U0| + |2;3;)(3;0|, 






dl- 


|3;0)(00|-|lA)(ltO|- 


|2t3i)(2r0h 


3T3;)(3tO| 


- 


m3;)(l;0|-|2;3;)(2;0|, 




(5) 



where \iaja') = |v) ® \ja') is a product state of single elec- 
tron states in which dots i and j are occupied by an electron 
of spin (T, a' =tj i- The state |0) corresponds to the case of no 
electron occupying any of the quantum dots. The minus signs 
in the expressions arise from the anti-symmetric property of 
the fermion system with respect to the exchange of two elec- 
trons. 

We may then decompose the Hamiltonian Hd into one- and 
two-electron parts 



H, 



Y^Ha+Y^Haa+Y.^- 



(6) 



a^a' 



where H^r represents the energy of the quantum dot system 
containing only a single electron of spin a, the Hamilto- 
nian Haa represents the case when two parallel-spin electrons 
are injected to the quantum dots, and Hera' (c 7^ cr') repre- 
sents the case of two opposite-spin electrons present in the 
quantum dot system. 

The single- and two-electron Hamiltonians can be conve- 
niently written in matrix forms. The space of the quantum dot 
system containing a single electron is spanned by three state 
vectors, |lcr), 12^), |3cr), the space of the system containing 
two parallel-spin electrons is spanned by three state vectors, 
|lcr2o.), |lcr3cr), |2cr3cr), and the space of the system contain- 
ing two opposite-spin electrons is spanned by nine state vec- 
tors, ll^U), |1^2J, |2tU>, |2t2;), 13^3;), |lt3;), 12^3^), 
|3fl4.), |3f24^). It is then straightforward to show that the 
Hamiltonian Ha written in the single-spin basis is of the form 



Ha 









Ti 





Al 


T2 


Ti 


T2 


A2 



Ell, 



(7) 



where Ai = E2 — Ei and A2 = E3 — Ei are energy differ- 
ences between the dot 1 and the dots 2 and 3, respectively, and 
a single energy factor Ei I, in which / is the unit matrix, has 
been pulled out for the convenience of diagonalization of the 
Hamiltonian and interpretation. This term can be referred as 
a reference energy about which the energies of the states are 
centered. 

The Hamiltonian Hera written in the basis of two parallel- 
spin electrons is of the form 



Ti -T2 

Ha.= \ Ti Ai+<5i 

-T2 A2 + S2 



-{Ei+E2 + Ui2)I, 



(8) 



where 61 = U13 - U12 and 62 = U23 - t^i2- 

Finally, the Hamiltonian Haa' written in the basis of two 
opposite-spin electrons is represented by the following 9x9 
matrix 
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(9) 



where 



(5ii = -Ai + L/11 - C/12, (522= Ai+[/22-f/l2, 

fe = A2 + A3 + f/33 - [/12, (5i3 = A2 + (5i, 

,523 = A3+J2, (10) 

with A3 = £^3 — £^2- 



C. Master equation 

We now suppose that the quantum dots and the NMR are 
weakly coupled to the electrodes and to the thermal reservoir, 
respectively. This "reduced" quantum system, i.e. the quan- 
tum dots plus the NMR can be described by a reduced den- 
sity operator p, which is obtained by tracing the density ma- 
trix of the entire system pT over the space of the Cjka and ciq 
modes fi23i. J27i - t29il 



= Tr 



(e,r)PT- 



(11) 



We shall derive the master equation for this reduced density 
operator by assuming that the reservoir modes and the elec- 
trons modes in the electrodes are (5-correlated, i.e. they appear 
as white (frequency independent) noises to the NMR and to 
the quantum dots, respectively. If we limit the interaction of 
the reduced system to the second-order in the coupling con- 
stants Vjk and cjq, we then find that under the Born-Markov 



and the rotating-wave approximations, the the reduced den- 
sity matrix p satisfies the master equation OOll 



dt 



p = -i[Hd + H„i,p\ + Cpp + CeP, (12) 



in which 



Cpp = x(n-p + l)7p {2dpa* — a*ap — pd*a) 

+ —np"ip {2d* pa — dd* p — pdd') , 
1 ^ 

i — l (7 
+ 2 H ^3'^ (24<tp4<t - dl^d'iaP - Pdlad-3aj , (13) 

(7 

where Cpp is an operator representing processes that lead to 
thermalization of the NMR phonon mode with the rate jp, 
and rip = [cxp{ojm/kBTp) — 1]^^ is the number of ther- 
mal phonons of frequency uj„i at temperature Tp. The opera- 
tor CeP represents the one-way tunnelling process of injection 
of electrons to the quantum dots 1 and 2 from the source elec- 
trodes with chemical potential /ii.2, and the damping of the 
dot 3 by the one-way electron tunnelling out to the drain elec- 
trode with chemical potential p^. The electrons tunnel into 
the dots 1 and 2 with rates Fi^ and F2ct, respectively, whereas 
they tunnel away from the dot 3 with the rate F3cr, as illus- 
trated in Fig.[T] 

The difference between the source and drain electrodes in 
the master equation ( fT2] i results from the difference between 
the corresponding chemical potentials, although the Hamil- 
tonian described in Eq. (|2]i is identical for all electrodes. In 
the derivation, we have assumed that the chemical potentials 
/ii.2 and /i3 obey the relations, /ii.2 > i?i,2 and p^ < E^, 
which are required for the tunnelling processes to occur in 
one direction. For these chemical potentials in the equilib- 
rium reservoirs, the Fermi distributions in the source and drain 
electrodes can well be approximated by /i(£'i) = 1 {i = 1,2), 
and fslE^) = 0. Thus, in the infinite-bias limit of /ii_2 — > oo 
and /i3 — > — oo, the three tunnelling rates have the form 



T,a^2nY,\V^k\^S{E,^e,k), 



(14) 



which is assumed to be independent of the energy Ei and 
spin a. 

D. Rate equation for the average phonon number 

Having derived the master equation for the NMR coupled 
to the quantum dots system, we now turn our attention to the 
derivation of a rate equation for the average phonon number 
in the resonator modes. In doing this, we will assume that the 
coupling between the quantum dots and the resonator is weak, 
cti.2 "C a;„i, Ti2- In such a case, we may use the second-order 
perturbation theory to eliminate the degrees of the freedom of 



the quantum dot system. A similar procedure has been used in 
the framework of single-electron dark states lUOl Ull Hsl Bill . 
We may extract from the Hamiltonian (|3) the part describ- 
ing the interaction between the dots 2 and 3 and the NMR 

Ha = a^{a + a^) {h2a +n3cr), (15) 



where, for simplicity, we put ai = 02 = a, and write the 
master equation (fT2l i in the form 



—p^ Cap + CaP, 



(16) 



in which the Liouvillian 

Cop = -i[Hd,p] - iuirn[a^ a, p] + CeP + Cpp (17) 

represents the dynamics of the quantum dots and the phonon 
mode in the absence of the electron-phonon coupling, and 



CaP = -i[Ha,p] 



(18) 



is the interaction between the quantum dots and the phonon 
mode. In the limit of a <C cdm, the interaction (fTSl l can be 
treated as a weak perturber to the system. 

We now follow the procedure of Cirac et al. 0211 of using 
the second-order perturbation theory in respect to the coupling 
strength a, to obtain from the master equation (fT&t a set of 
rate equations for the occupation probabilities p„ = {n\p\n) 
of the NMR phonon state with n excitations. Straightforward 
but lengthly calculations yield the following equation 

-TfPn = [{flp + l)7p + ^-] [{n + l)P«+i - npn] 

+ (riplp + A+)[np„_i - (n + l)p„], (19) 

in which the transition rates A± are given by 

A± = 2a^RcS{±uj„,) 



2a^Rc / dT{d{T)d{0))sC 



^lUJmT 



(20) 



where {O{t)O{0))s is the correlation function for the quan- 
tum dots operators 0{t) = J^A^^^'i''') + "-So- (''')]' evalu- 
ated in the absence of the coupling a. The two-time correla- 
tion function is evaluated using the quantum regression theo- 
rem 1I33I1 . which states that the two-time correlation function 
{0{t)0{Q))s satisfies the same equations of motion for t > 
as the corresponding one-time averages {0{t)). In turn, the 
one-time averages can be expressed in terms of the matrix el- 
ements of the reduced density operator of the the quantum 
dots pd = Tr,.p, where the trace is taken over the space of the 
NMR mode. 

To examine the occurrence of cooling we look at the aver- 
age number of phonons in the steady-state of the NMR mode, 
(?i) — J2n "Pn- Using Eq. ( fT9] l, we obtain the following sim- 
ple equation of motion 



dt 



(n) = - (7p + ^- - A+) (n) + -fpfip + A^ 



(21) 



In order to calculate (n) we need to determine the rates A± . It 
is seen from Eq. (ISTT i that A^^ and A^ are heating and cooling 
rates, respectively. They represent processes which can in- 
crease (heating) and decrease (cooling) of an excitation of the 
NMR mode when at most two electrons are simultaneously in- 
jected to the quantum dots. Under the condition of A _ > A+, 
the steady-state average number of phonons is of the form 



{n)s = 



IpUp ■ 



■A4 



7p- 



A^-A, 



(22) 



It is clear by examination of Eq. (l22l i that a significant cooling 
of the NMR can be achieved when A^ ^ 7^, A^. 

In the following section, we shall use Eq. (l22l) to study the 
conditions for cooling of the NMR to its ground state with the 
help of the two-particle dark state. As we shall see, electrons 
trapped in the two-particle dark state can be used as the carri- 
ers to transfer the phonon energy from the NMR mode to the 
drain, thereby cooling the NMR to its ground state. 



III. DARK-STATE COOLING OF THE NMR 



Now we demonstrate how to apply the two-particle dark to 
cool the NMR to its ground state. The task is to determine the 
mean phonon number in the steady-state which expresses the 
information about temperature of the oscillating resonator To 
see the advantages of using the two-particle dark state, we first 
briefly consider the cooling mechanism with a single-particle 
dark state Il20ll2ill3ll35ll. 



A. Cooling with the single-electron dark state 

In the strong Coulomb-blockade regime, at most one elec- 
tron can tunnel into the quantum dot system. With only one 
electron present, the quantum-dot system is described by the 
Hamiltonian He alone, Eq. (|7]i. We may diagonalize the 3x3 
matrix appearing in Eq. (|7) to find the effective (dressed) 
states of the electron. On calculating the eigenvalues and 
eigenvectors of the matrix one can find that all of the eigen- 
values are different from zero and the corresponding eigenvec- 
tors all involve the state |3ct). This indicates that in general no 
electron trapping (dark) state is created because the electron 
can escape from the dot 3 to the drain electrode. However, it 
is easily verified that in the limit of Ai = 0, the eigenvalues 
of the matrix are 



1 



Ai = 0, A± = ^A2±W(r2 + r|) + iAi, 



and the corresponding eigenvectors are of the form 



(23) 



\^ia) = cos?7|U) - sin?/|2^), 

\^+a) = cos6'|3<j) +sin6'(sin?7|lcr) + cos?7|2o.)), 

\ip^a) = sin6'|3cr) -cos6'(sin?7|lo.) +cos?7|2o-)), (24) 



where tan 77 = T1/T2, and 



tan6 



Q 



n 



(25) 



with n = ^Al + 4(r2 + Tl). 

We see that one of the eigenvalues (Ai) is zero and 
the eigenvector corresponding to Ai does not involve the 
state |3cr)- This shows that by choosing two pumped dots of 
the same energies, E2 = Ei, we may achieve a dark state 
that is decoupled from the dot 3. This statement remains true 
even if A2 = 0. Note that the creation of the dark state is 
independent of the ratio T1/T2. 

It is straightforward to see from Eq. (l2Jt that A+ > Ai > 
A_. Thus, the state \f+cr) is shifted upwards from the dark 
state energy by the amount A+, whereas the state \^p~cr) is 
shifted downwards from |(/?icr) by the amount — A_. Hence, 
the electron will experience gain of its energy when making a 
transition from the dark state \Lpicr) to the state \f+a) and ab- 
sorption of the energy when making a transition from the dark 
state to the state |(^-(t). Thus, cooling of the NMR is possible 
only when the electron interacting with the NMR makes tran- 
sition to the upper state. When A2 = 0, i.e. all three dots are 
identical, the states |i^+c,) and \'P-a) are then symmetrically 
located about the dark state |<picr). In this case, the transi- 
tion frequencies |i^icr) -^ W+a) and \tfia) — > \f-<y) over- 
lap. Consequently, the cooling and heating frequencies over- 
lap resulting in no effective reduction of the average number 
of phonons. Hence, the one-electron cooling can be observed 
only when A2 7^ 0, i.e. when the states \f+a) and \f~cr) are 
not equally shifted from the dark state \ificr)- 



mK, corresponding to fip = 21, equal couplings Ti = T2 = T 
and several different values of the detuning A2. The mini- 
mum of {n)s is seen to occur at the frequency cUm = A2/2 + 
y^2T^ + A2/4, which corresponds to the NMR frequency on 
resonance with the transition from \ipia-) to the state \f+a)- In 
the case of large energy difference between dots 3 and 1, the 
mean phonon number {n)s can be reduced significantly. For 
example, when A2 = 2r, the minimum value of the mean 
phonon number can reach 0.006, so that we may speak of 
cooling the NMR to its ground state. Figure |2] also shows 
that a significant reduction of {n)s is achievable in principle 
at large A2. When A2 = 0, no significant reduction of (ri)^ is 
observed which makes us to conclude that no cooling can be 
achieved with three identical quantum dots. 





14.5 15.0 15.5 16.0 16.5 17.0 



FIG. 2: The steady-state average number of phonons {n)s plotted as 
a function of the mechanical frequency ojm for a — 2r, an initial 
temperature Tp = 100 mK (corresponding to rip — 21), 7^ = 2 x 

lo^^r, Ti = r2 = r, Ti = T2 = lor with r^ = o.sr and 

different A2: A2 = O.SF (solid line), A2 = T (dashed line), A2 = 
2r (dotted-dashed line) and A2 = 4r (dotted line). 

The above considerations are illustrated in Fig. |2] which 
shows the steady-state average number of phonons as a func- 
tion of the frequency oj™ for an initial temperature Tp — 100 



no. 3: The mean phonon number {n)s as a function the mechanical 
frequency ujm for a = 2r, an initial temperature Tp — 100 mK 
(corresponding to rip = 21), Fi = r2 = P, 7p = 2 x 10~*F, F3 = 
0.5r, A2 = F, the average tunnelling amplitude T = (Ti+r2)/2 = 
lOF and different ratios /3 = T1/T2: ^ ^ I (solid line), /3 = 4 
(dotted line), /3 — 9 (dashed line) and 13 — 19 (dotted-dashed line). 



It is interesting to extend the present analysis to the case of 
tan 7] y^ 1, corresponding to a dark state with unequal ampli- 
tudes of the superposed states |lcr) and |2cr). A brief look at 
Fig-ffland Eq. (l24l i suggests that a dark state with tan 77 > 1 
would be more suitable for cooling since the asymmetry in the 
population of the dots 1 and 2 with P22 > pii could lead to a 
stronger interaction between the dot 2 and the NMR. Unequal 
tunnelling rates with Ti > T2 is the mechanism for creat- 
ing P22 > pii- This is shown in Fig. [3] where we plot the 
steady-state mean phonon number {ji)s as a function of curn 
for different values of the ratio /3 ~ Ti/T2. It is seen that the 
effect of going from /? = 1 to /3 > 1 is to reduce rather than to 
increase the cooling efficiency thus leading to the growing up 
of the minimum of phonon number The reason is that unequal 
tunnelling rates are equally effective in destroying the asym- 
metry between the cooling and heating frequencies. These 
considerations imply that the best conditions for cooling with 
the single-electron dark state are Ti — T2 = T and A2 ^ F. 



B. Cooling with the two-electron dark state 



which in the symmetric case simplifies to 



We now consider the cooling mechanism when two elec- 
trons can be simultaneously injected into the quantum-dot sys- 
tem. Let us first specialise to the case when only two parallel- 
spin electrons are present. This situation is described by the 
Hamiltonian ^, which in the case of Ei = E2 = E and 
5i = 62 ^ S can be written as 



H„ 






Ti -T2 ^ 




Ti 


A2+6 


+ {2E+Ui2)I 


-T2 


A2+d^ J 





(26) 



The 3x3 matrix appearing in Eq. (|26> can be diagonalised and 
the dressed states of the electrons can be identified. However, 
one can easily verified, without going into the detailed deriva- 
tion, that the dressed states all would involve the state |3cr), 
which is dissipative due to its coupling to the drain. This dis- 
sipative state is involved regardless of whether the dots are 
identical (A2 = 0) or not (A2 7^ 0). Thus, no trapping of 
the electrons could be achieved. Consequently, the electrons 
would escape to the drain without interacting with the NMR. 
As a result, no cooling of the NMR could be observed in the 
interaction with two electrons of the same spin. 

We now generalize the situation to the case of two electrons 
of opposite spins. To see the advantage of the two-electron 
over the single-electron case, we assume that dots are identi- 
cal, i.e. Ai = 0{i = 1, 2, 3). Notice, that no single-electron 
cooling was predicted in this regime. Therefore, we may dis- 
tinguish between one- and two-electron coolings. The case 
of two electrons of opposite spins is described by the Hamil- 
tonian ^ which, in general, has a complicated form but it 
simplifies in various special cases. We consider two of these 
cases. A symmetric case in which we choose the charging po- 
tential energies Uu = U22 = t^33 and t/13 = C/23 = t/12, and 
an asymmetric case of C/n = U22 7^ C^33 and [/13 = C/23 ¥" 
Ui2- In the symmetric case the de tunings satisfy the condition 
<^ii = <^22 = 1^33 = Su and ^13 = i523 = 0, whereas in the an- 
tisymmetric case (5ii = (522 = <5„, and S13 — S23 — S^j 7^ 0. 
We shall assume additionally that Ti = T2 = T. In each case 
we will be able to find analytical formulae for the two-electron 
dressed states of the coupled quantum dots. 



/i5„ 





H„ 



T T \ 
T T 
T T 
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{2E + Ui2)I- 



(27) 



When the detuning (5u = 0, a diagonalization of the ma- 
trix (IZTT i leads to dressed states | $) which are related to that 
in the bare-state basis by the transformation 



|*)=W,|*) 



with the corresponding eigenvalues 



(28) 



Ai = -2V2T, X2 = 2V2T, A3 = A4 = -V2T, 
As = Ae = V2T, A7 = Ag = Ag = 0, (29) 

where \^) is a column vector composed of the nine bare states 



I*) = [|ltU), |lt2;), |2^U>, |2^2^), |3t3;), |lt3;) 



2t3_L}, 3tl 



^f^i 



t^i/ 



3t2 



tH> 



(30) 



and Ws is a 9 X 9 transformation matrix whose explicit form 
is given in the Appendix lAl 

A careful examination of the transformation matrix W^, 
Eq. (lAlb . reveals that among the dressed states there is a 
state |$g) which is a linear superposition involving only bare 
states of the dots 1 and 2: 

\^9) = I (|ltU> - |lt2i) - |2tli> + |2t2i» . (31) 

All the remaining states involve bare states of the dot 3. Since 
the electrons can tunnel to the drain electrode solely from the 
dot 3, this implies that the state j^g) is an example of a two- 
electron trapping (dark) state. In other words, two electrons, 
once trapped in the state |$g), cannot tunnel into the drain 
electrode. 

We now couple the resulting electron dress states to the me- 
chanical resonator to see if the presence of the two-electron 
dark state dSTT i could lead to cooling of the resonator. The 
interaction between the quantum-dot system and the NMR is 
determined by the Hamiltonian (fTsT l. We make the unitary 
transformation 



1. The symmetric case 



The contribution of two opposite-spin electrons to the dy- 
namics of the system is described by the Hamiltonian (|9|l. 



Ha{t) ~ cxp{iHot)Ha cxp{—iHot), 



where 



i^o 



OJr, 



H, 



(32) 



(33) 



and find that the Hamiltonian ( fTSl l in the dressed-state basis 
takes the form 



t{'>) 



H':>it) = 



_ " r.t 



a' 



*5) + |*6))($9|e-'("™-^^'* + •••] + H.c.} 



(34) 

In writing the expression (|34] | we have limited ourselves to 
those terms that correspond to transitions from the dark state 
|$g) and we only concern the case that the frequency of the 
NMR equals to \/2T. The remaining terms corresponding to 
transitions occurring at different frequencies can be neglected 
for their rapid oscillation. Another kind of terms involving 
dressed states which contain the state |3cr) can be neglected 
because the populations on these states are negligible. 




FIG. 4: Schematic diagram of possible cooling and heating transi- 
tions from the dark state |"1>9) to other dressed states of the system 
in the case of the symmetric coupling between the dots. The absorp- 
tion (cooling) processes correspond to the transitions indicated by 
solid arrows while the emission (heating) transitions are indicated by 
dashed an'ows. 

The Hamiltonian ( |34] | determines which of the transitions 
can occur with the absorption of a phonon from the NMR 
mode and which with the emission of a phonon to the NMR 
mode. Figure |4] shows the two-electron dressed states of the 
quantum dot system with possible transitions from the dark 
state I $9). The electrons trapped in the state |$g) could make 
transitions with frequency aj,„ = \/2T to the upper states |<I>,5) 
and I $6) by absorbing a phonon from the NMR mode. Of 
course, the absorption of the phonon would result in a cooling 
of the NMR to lower temperatures. Unfortunately, the trapped 
electrons could also make transitions with the same amplitude 
and at the same frequency to the lower states [$3) and [$4) by 
emitting a phonon to the NMR mode. Since both processes, 
the absorption and emission of phonons occur with the same 
amplitudes and at the same frequencies, no net cooling of the 
NMR could be achieved. In the following, we will look into 
the cooling mechanism subject of an asymmetry in the charg- 
ing potential energies coupling the injected electrons. 



2. The antisymmetric case 

Our interest is to achieve cooling of the NMR with the help 
of the two-electron dark state. In order to study this problem, 
we try to unbalance the symmetry between the absorption and 
emission processes from the dark state. This could be done, 
for example, by introducing an asymmetry between the charg- 
ing potential energies. 

Suppose that Un = U22 ^ U33 and U13 = U23 ^ Uri- 
In this case, (5i3 = (^23 = Sw 7^ 0. We then find that the 
Hamiltonian (|9]i reduces to 

/(5„000 0r0T0\ 
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(35) 

A diagonalization of the Hamiltonian ( l35T l produces com- 
plicated formulae for the eigenvalues (energies) and the eigen- 
vectors (dressed states) of the system. Even at (5„ = there 
are still three parameters involved in these formulae, T, 833 
and Sw Therefore, for a simplification of the formulae, we 
express 833 and 5^ in terms of T. For example, setting the 
detunings J33 = 5u + 2T, 6^ = T and 6u = 0, we find simple 
expressions for the eigenvalues 

Ai = AT, A2 = A3 = 2T, A4 = As = T, 

\e = -2T, Xr = 0, Ag = Ag = -T, (36) 



and the corresponding dressed states 

|*)=WJ*), 



(37) 



where Wq is a 9 x 9 transformation matrix whose explicit 
form is given in the Appendix lAl 

An inspection of Eq. (IA2I) reveals that one of the dressed 
states, I $7), is of the same form as the state ( ISTT i. i.e. the 
state |$7) is a dark state involving bare states ofonly the dots 1 
and 2. We have checked that independent of the choice of S33 
and 6w, there is always an asymmetric dressed state that in- 
volves bare states of the dots 1 and 2 only. Hence (5„ = is 
the general condition for the creation of the dark state in the 
system. 

Given the electron dressed states, we may transform the in- 
teraction Hamiltonian Ha into the dressed-state basis and find 



m 



r)(t) = ^{at[(|ci>g) + |cf>,))(cf>,|c' 



{uj^-T}t 



R.c}, (38) 



a 

'V2 

where as before for the symmetric case, we have explicitly 
listed only the dominant terms in the coupling of the dressed 



system to the NMR. One can notice that these dominant terms 
involve operators describing transitions to and from the dark 

state I $7). 




\7.T 



FIG. 5: Two-electron dressed states of the quantum-dot system in 
the case of the asymmetric coupling between the dots. The solid ar- 
rows indicate dominant cooling transitions from the dark state |<1>7) 
to the upper states |$2) and l^a), and dashed arrows indicate emis- 
sion (heating) transitions to the lower energy states l^g) and |<1?9). 

Figure |5] shows the two-electron dressed states of the 
quantum-dot system for an asymmetric coupling between the 
dots. The arrows indicate transitions from the dark state that 
are induced by the coupling of the dressed system to the NMR. 
We see that the dark state |$7) is strongly coupled to only four 
states of the dressed system, two upper and two lower energy 
states. Thus, similar to the symmetric case both cooling and 
heating processes can occur However, a notable distinction 
between the symmetric and antisymmetric cases is that in the 
asymmetric case the absorption and emission processes oc- 
cur at different frequencies. It is easily verified from Eq. (l38l l 
that the absorption processes occur at frequency ujrn = ST, 
whereas emission processes occur at ujm = T. This means 
that at the frequency a;„i = 2T, the absorption of a phonon 
from the NMR is not accompanied by an emission of a phonon 
of the same frequency. As a consequence the NMR can be 
cooled to its ground state. 

This is illustrated in Fig. |6] which shows the steady-state 
mean phonon number (77,) ^ as a function of the frequency a;,„ 
of the NMR for different tunnelling rates of the electrons to 
the drain. Cooling of the resonator shows up as a reduction 
of the mean phonon number with a clear minimum at Wm = 
2r. The minimum of (n)s, corresponding to the maximum 
of the cooling of the NMR, is achieved at small values of the 
tunnelling rate Fa and the effect of increasing F3 is seen to 
reduce the cooling efficiency. Taking into account the mean 
phonon number of rip ~ 21 in the initial equilibrium thermal 
state, the value of (71)3 ~ 1.3 x 10^'^ achieved with the two 
electron dark state means that the NMR has been cooled to 
its ground state. With the practical values of the frequency 
Wm = 100 MHz, the predicted minimum value of (n)^ ~ 
1.3 X 10"'^ corresponds to a temperature of the NMR reduced 



to T„, ~ 0.72 mK. 
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FIG. 6; Dependence of the steady-state mean phonon number (n)s 
on frequency ojm for the asymmetric case of 5™ = T with an initial 
temperature Tp = 100 mK (corresponding to n^ = 21), Fi = r2 = 
F, r = 10F,a = 2F, (5„ = 0, A^ = (i = l,2,3),7p = 2 x 
10~''r,(533 = 2r, and different tunnelling rates F3: F3 = O.SF 
(solid line), r3 = F (dashed line), F3 = 2r (dotted-dashed line). 

Figure|7]shows the effect of the tunnelling rate Fs, at which 
electrons escape from the system, on the magnitude of the 
minimum value of (n)s. The effect of increasing Fs is clearly 
to reduce the minimum value of the mean number of phonons. 
For small tunnelling rate T the increase of the minimum value 
of (n)s with Fa is considerably more rapid than it is for 
large T . This is readily understood if one recalls that increas- 
ing of the escape rate of the electrons results in a shortening 
of the interaction time of the electrons with the NMR. The 
shorter interaction time implies less of the energy taken from 
the NMR. Thus, the minimum value of (n)^ degradates with 
the increasing tunnelling rate F3. 

The above analysis of the cooling mechanism were based 
on the energy structure of the system of the coupled dots. Fur- 
ther insight into the role of the two electron dark state in the 
cooling process is gained by considering the dynamics of the 
dressed states 1$;). Because it is precisely the effect of trap- 
ping of the population in the state [$7) that was crucial for the 
cooling mechanism, we derive the equation of motion for the 
population of the state [$7) that is represented by the density 
matrix element ,077. To illustrate that the dressed states that 
involve states of the dot 3 cannot be trapping states, we also 
derive equations of motion for the dressed states |<1>2) and [$3) 
that according to Eq. ( |37] | are linear superpositions involving 
the state |3cr). 

Consider the evolution of the quantum dot system alone, i.e. 
in the absence of the NMR. The evolution is determined by the 
master equation that is obtained from Eq. (fT2l i by tracing over 
the NMR states 



d 
JtP 



-i[Hd,p] + CeP, 



(39) 



where Hd is given in Eq. ^ and £,./? is given in Eq. ( fT3] ). 
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FIG. 7: The dependence of minimum value of {n)s on the tunnelling 
rate Fa for the asymmetric case of 5^ — T with an initial temper- 
ature Tp = 100 mK (corresponding to rip — 21), Fi = r2 = F, 

= 2r, Q = 2F, 5„ = 0, Ai = (i = 1, 2, 3), jp = 2x lO^^F, 



^33 = 2T, and different T: T = 5F (solid line), T - 
line), T = 15F (dashed-dotted line). 



lOF (dashed- 



We conclude this section by pointing out that the condi- 
tions predicted by Fig. |6] for cooling of a NMR with the two- 
electron dark state may be met in the current experimental sit- 
uations. The value of the tunnelling rate r/27r = 5 MHz used 
to generate Fig. |6] corresponds to the typical values occurring 
in electron transport experiments where the tunnelling rates 
ranging from 10 kHz to 10 GHz were measured ll36ll37ll . The 
predicted cooling of the NMR occurs at frequencies Wm/F w 
20 corresponding to the value of ujm ~ 100 MHz, which ap- 
pears to be practical OSll . Note further that the expected cool- 
ing times tc correspond to that required for the mean phonon 
number to reach the steady-state, which according to Eq. (ISTT l 
are of the order of (7p + A_ — A+)~^. With the parameter 
values of 7j, = 2 X 10~^r, corresponding to a realistic quality 
factor Q - 10^ iHliSl, ^- - 4.577r and A+~5x IQ-^T, 
the cooling could be achieved at times tc — 4.4 x 10"^ s. 

Finally, interesting further generalizations of the results 
presented in this paper could include other configurations of 
the quantum dot system such us l/-type or linear systems. It is 
also worth considering systems composed of a large number 
of regularly distributed quantum dots. 



When we project the master equation (1391 1 onto \^i) (i = 
2,3,7) on the right and ($i | on the left, we obtain the follow- 
ing equations of motion 

P77 == r (p<^i^ij, + Pipi^i J - -^iSu {pn + 2/957 + 2/967) , 

2 1 

P22 = -0^3/322 + 7I 

3 6 



(x/2-l) p^_^_ 



2 1 

P33 = -XI3P33 + tI 

3 6 



2P¥'uu + (1 + \/2J p^^^^, 

i:i&u P23 + V2(p28 + P29) 
2P¥'im + (1 + \/2J p^^j^^^ 



(x/2-1 



Pv-i-^ 



-i^^. 



P23 + V2(p38 + P39) , 

(40) 



where Pifi^i^ is the population of the single electron 
state |(/3icr) and py = Imp.y is the imaginary part of the co- 
herence between states i and j. 

We see from Eq. ( |40] | that the dressed state [$7) does not 
decay, but can be populated by transitions from the single 
electron states |(^i^) and \'fi{). In contrast, the states [$2) 
and I $3) that involve the state of the dot 3 decay with the 
rate 2r3/3. The population state p-ji is coupled to the other 
state through the detuning (5„. In the steady state {pu = 0) 
with 6u = 0, the population can be completely transferred 
to the state |$7). Thus, the equations of motion (|40] | clearly 
demonstrate that in the steady-state with 5u = 0, the state | $7) 
becomes a trapping state. 

The expressions ( l40t also show that the single electron 
states play the major role in the two electron cooling process. 
The dark state |<I>7) is populated with the rate F from the single 
electron states |</Jif ) and \fii)- When 5u ~ 0, no population 
can be transferred from the other two-electron states. 



IV. SUMMARY 

We have studied the dark state cooling technique in a sys- 
tem composed of three coupled quantum dots interacting with 
an NMR. Two cases of coupling of the NMR to single and 
two-electron dark states have been considered. Cooling of 
the NMR to its ground state can occur in both cases. In the 
case of the interaction of the resonator with a single-electron 
dark state, no cooling of the resonator occurs unless the quan- 
tum dots are not identical. The steady-state cooling is possible 
only if the energy state of the quantum dot coupled to the drain 
electrode is detuned from the energy states of the dots coupled 
to the electron source electrode. In the case of identical dots 
the cooling and heating processes were found to occur with 
the same amplitudes and at the same frequency. This results 
in no effective cooling of the NMR. When the dots are not 
identical, the phonon emission (heating) and phonon absorp- 
tion (cooling) transitions occur at different frequencies giving 
rise to an effective cooling of the NMR. 

In the case of the interaction of the NMR with the two- 
electron dark state, we have showed that an effective cooling 
can be achieved even with identical quantum dots subject of 
asymmetries only in the charging potential energies coupling 
the injected electrons. In the symmetric case, we have shown 
that the cooling and heating transitions occur with the same 
amplitudes and at the same frequencies. An asymmetry in the 
charging potential energies is found to lead to unequal shifts 
of the effective dressed states of the systems which results in 
the cooling and heating frequencies to become detuned from 
each other. This shows that the physics of the two-electron 
dark state cooling process is quite similar to that of the single- 
electron dark state cooling. However, there is an important 
difference that these two cooling processes occur at signifi- 
cantly different frequencies. This indicates that the frequency 
at which the resonator could be cooled to its ground state can 
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be changed by switching from the one-electron to the two- 
electron Coulomb blockade process. 
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In this Appendix we give the transformation matrices from 
the bare electron states to the two-electron dressed states. In 
the case of the symmetric coupling the transformation matrix 
is of the form 



W, 
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whereas in the case of the asymmetric coupling, the matrix is given by 
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